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1. INTRODUCTION 

Present and future lattice calculations involv- 
ing b and c quarks include some of the most im- 
portant applications of lattice gauge theory to 
standard model physics. These include the heavy 
meson decay constants, the BB mixing ampli- 
tude, and various semileptonic decay amplitudes, 
which are all crucial in extracting CKM angles 
from experimental data. They also include the 
extraction of a s from the charmonium and bot- 
tomonium spectra. 

Bound states of heavy quarks and antiquarks 
(quarkonia) have another crucial role to play in 
the development of lattice gauge theory: they 
provide systems in which the estimation of the 
errors inherent in current lattice calculations can 
be done in a more reliable and robust way than 
is possible for the light hadrons. The reason is 
that the quarks in these systems are relatively 
nonrclativistic. Coulomb gauge wave functions 
calculated on the lattice may be used to aid in 
the estimation of finite volume and finite lattice 
spacing errors, and of the effects of quenching. 
We have a much better idea of what to expect 
in lattice calculations of these systems since po- 
tential models may be used to obtain the leading 
behavior in i> 2 /c 2 .jl| 

Chris Sachrajda and I will split the subject of 
heavy quark physics in these proceedings. His re- 
view H will concentrate on the part of the subject 
which involves the weak interactions. Mine will 
concentrate on the part which does not. 



2. LATTICE FORMULATIONS 
HEAVY QUARKS 
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When the lattice spacing a is smaller than the 
Compton wave length of the quark 1 / m, the stan- 
dard relativistic action of Wilson may be used. 
Cutoff effects may be removed by taking the cut- 
off l/a to infinity. The bare lattice action may 
also be viewed as an effective field theory of QCD 
at the cutoff scale. Cutoff effects in an effective 
field theory are removed by adding higher dimen- 
sion interactions to the bare Lagrangian while 
keeping the cutoff fixed. To remove the effects 
of the cutoff to a finite order in a, a finite number 
of interactions may be added to the bare lattice 
Lagrangian. When oAqcd "C 1, perturbation 
theory may be used to calculate the required co- 
efficients of the new operators. Q The ability to 
remove cutoff effects pcrturbatively will probably 
be spoiled eventually, perhaps at a small power of 
a due to effects presently not understood, almost 
certainly at a relatively large power of a due to 
instantons. 

The dynamical scales in bound states are small 
compared to the fermion mass in QED and in 
QCD for the c and b quarks. It is often advan- 
tageous in these systems to formulate the field 
theory nonrelativistically as an expansion in l/m 
|] ||], keeping the cutoff at or below m. The non- 
leading terms in nonrelativistic expansions have 
dimension higher than four. Loop corrections in 
these effective field theories diverge if the cut- 
off is removed. Cutoff effects must be removed 
by adding higher dimension interactions to the 
Lagrangian or by raising the cutoff to the new 
physics scale (m), switching to the relativistic, 
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renormalizable version of the theory, and then 
taking the cutoff to infinity. 

When the kinetic energy of the heavy quark is 
small compared to the typical interaction energies 
(as it is in bound states containing a single heavy 
quark), the kinetic energy may be treated as a 
perturbation. In this static approximation 0j6|, 
the lowest order fermion action is just 



^static = 4>*iDt 



(1) 



and the unperturbed quark propagator is just the 
timclike Wilson line. In the general case, includ- 
ing quarkonia, the lowest order potential and ki- 
netic terms of the Lagrangian of Nonrelativistic 
QCD (NRQCD) §,§ (the terms on the first line 
of Equation ||) must be included in the unper- 
turbed Lagrangian. 



'{iD t 



+ 2m 

+ £*-B } , 



(2) 



Higher order terms in may be added as per- 
turbations. 

In processes such as the semileptonic decay of 
heavy-light mesons, in which one heavy quark de- 
cays into another lighter (but still heavy) quark 
with a high velocity relative to the first, it is pos- 
sible, and useful, to formulate the static approxi- 
mation and nonrelativistic QCD as expansions in 
the small internal quark momentum around some 
large, external meson momentum. [H Lattice im- 
plementations of this idea have been proposed |n| 
and are reviewed by Sachrajda.|| 

2.1. The improvement program for Non- 
relativistic QCD 

In a recent paper 11 1, Lepage et al. have sys- 
tematically examined the improvement program 
for NRQCD with the goal of reducing systematic 
sources of error from all sources to under 10%. 
This program involved the following elements: 

1) Since NRQCD has been formulated as a non- 
renormalizable effective field theory, cut-off ef- 
fects are removed not by taking the cut-off to in- 
finity, but by adding additional operators to the 
bare Lagrangian (the dots in Eq. 0) . The infinite 



number of possible operators must be ordered ac- 
cording to expected size of their effects on the 
physics. For heavy-light systems, the operator 
ordering is simply an expansion in 1/ra, that is, 
in the dimensions of the operators. For heavy- 
heavy systems like quarkonia, the expansion is 
complicated by the presence of large quark ve- 
locities which do not fall to zero with the quark 
mass. Operators with the same dimension (such 

as fm arL< ^ 2m a ' ^ ) are suppressed in their ef- 
fects on the physics by different powers of v (by 
v 2 and v 4 , respectively, in this case). 

2) Once the operators required for a given accu- 
racy have been established, their coefficients must 
be determined by requiring that the NRQCD La- 
grangian reproduce the Green's functions of ordi- 
nary QCD to this accuracy. 

3) Discrete forms of the required operators 
must then be defined. As with light quark ac- 
tions, finite a errors must be estimated. If nec- 
essary, correction operators Q must be added to 
the action. 

4) The coefficients of the operators arc modified 
by quantum effects. Many corrections have been 
calculated in mean field theory. [jll| The correc- 
tions for the quark energy shift, mass renormal- 
ization, and wave function renormalization have 
been calculated in full one-loop perturbation the- 
ory. Deviations between the mean field and 
one loop results are rather small, from 0-10%. 

The result is a systematic correction program 
in v, a, and a s . The correction operators in v 
and a may be included directly in the simulation 
action, or evaluated as perturbations using lattice 
or potential model wave functions. 

2.2. A New Action for Four Component 
Fermions 

Because coefficients of higher terms in the 
NRQCD Lagrangian such as are explicit func- 
tions of 1/m, the quantum corrections described 
in 4) above are also explicit functions of 1/m. 
These begin to diverge as ma is reduced below 
a value of order one, making the nonrelativistic 
expansion impractical. The Wilson action like- 
wise has been thought to have finite lattice spac- 
ing errors of order ma which blow up as ma is 
raised above one. Since the masses of the b and 
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c quarks are such that ma is often O(l) at cur- 
rent lattice spacings, calculations of such crucially 
interesting quantities as the heavy meson decay 
constants fg and fa have often involved awk- 
ward interpolations between results in the static 
approximation and results using Wilson fermions 
through a region where neither approximation is 
well behaved. || While such an approach is proba- 
bly workable, it would clearly be desirable to have 
a method for lattice fermions which did not begin 
to break down right in the region of interest. 

To approach such a method, we consider a lat- 
tice version of £ NRQCD with a few minor modifica- 
tions. Like Wilson fermions (tp), the fermions of 
NRQCD contain four components per site: a two- 
component quark field (</>) and a two-component 
antiquark field (\)- The bare mass is conven- 
tionally omitted in NRQCD calculations, but we 
are free to leave it in the theory. The usual Dirac 
coupling between quarks and antiquarks is absent 
(having been transformed into higher derivative 
interactions by the Foldy-Wouthhuysen transfor- 
mation), but we may add back a sufficiently sup- 
pressed amount of this interaction without spoil- 
ing the theory. We thus consider the following 
Lagrangian: 

C = 0*( Cl A-+m o -|^A+A-)0„ 

i 

+ C 3 <j>* ^ Vi^tXn 
i 

+ x*(-ci A++m -|^A+Ar) x „ 

i 

+ c 3 x* yVjA^, (3) 

i 

When ci = 1 (times a correction factor when 
ma ^ 1), C2 = — , and c 3 is negligible, it is a 
good Lagrangin for NRQCD. The point of writ- 
ing the NRQCD Lagrangian in this particular 
form is that the action becomes precisely the stan- 
dard Wilson action with the choice of parameters 
ci = C2 = c 3 = 1. It is thus possible to adjust the 
parameters in such a way that as mo is reduced, 
instead of blowing up, the theory turns smoothly 
into the Wilson theory. 

It is illuminating to expand the equation for 
Wilson propagators nonrelativistically when the 



mass is large. After normalizing the fields by 
^ 1 1 _ 6k fl^l (not -7= as is conventional) one may 
obtain 

6 nQ = [-£ + M + (1 - Ut fi ) (4) 

2 \m (1 + mo) (2 + mo) J ^ 

where £ is the energy eigenvalue obtained from 
the transfer matrix and M. = E v i = o = ln(l + mo). 
This is a lattice Schrodinger equation not un- 
like the one obtained from NRQCD, but it has 
some unusual features. Most important, the two 
"masses" in the equation, M. = Ln(l + mo) and 
tt = — + 7t~; tttt \ , are completely diffcr- 

M m a (1+mo) (2+m ) ' ^ J 

ent. If M. is used to fix the fermion mass when 
am 3> 1, the dynamically more important mass 
condition d£ / dp 2 = will be completely incor- 
rect. 

Kronfcld showed in his talk at this conference 
plf| that the two masses can be put back into 
agreement with the use of the action 

s = y^[-yi n y> ra 

n 

+ K t ^„(l - 7o)[/„,o^ I1+ 6 + h -c. 

+ k s ^J-0„(1 - 7i)J7„,iV n+ | +h.c] (5) 

i 

Thus, it seems that an action closely related to 
the Wilson action is a member of the class of 
actions suitable for NRQCD. One can go even 
further. In NRQCD and in the static approxi- 
mation, M. plays no dynamical role. It can be 
ignored, and is conventionally thrown away. This 
suggests that the standard Wilson action itself 
can be used when am > 1 as long as M. is ig- 
nored and d£/dp 2 — is used to fix the quark 
mass, as is done in NRQCD. 

This proposal is obviously correct in free field, 
where we can calculate the behavior of quark 
propagators exactly to see that the proposed in- 
terpretation makes sense. It is certainly correct 
in mean field theory, too. Mean field improve- 
ment of these fermions, as of Wilson fermions, is 
simply the absorption of a "mean link" uo (see 
Appendix |b|) into an effective k = uok and then 
proceeding as with free field theory. (A plausi- 
ble estimate of the mean link in this context is 
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probably uq ~ l/8/t c .) It remains to be shown 
whether the theory is somehow spoiled by renor- 
malization. 

Perturbatively, Green functions must be ex- 
panded in p 2 and a s ■ Each term in the expansion 
is an explicit function of the quark mass, since the 
theory must be solved exactly in ma. (The is also 
the case for the loop corrections of NRQCD. Jl| 
If these functions become singular or badly be- 
haved in some way, the theory could conceivably 
break down. The one loop perturbative correc- 
tions contain all of the ugliest features of Wilson 
and NRQCD perturbation theory simultaneously, 
and have only been begun. There is, however, one 
numerical calculation by El-Khadra indicat- 
ing that nothing too surprising occurs. The one- 
loop correction to the local current normalization 
for Wilson fermions with the naive normalization 
is M 



1 



(6) 



2k(1- 0.17c/ 2 )' 

The correct normalization with mean field im- 
provement is 



i 



(i-l^a-o.oe^) 



(7) 



The remaining perturbative correction, 0.06<7 2 , 
becomes an explicit (so far uncalculated) function 
of m (or k) in the new formalism which must not 
become singular if the theory is to make sense. 

Fig. ^] shows Eqs. || (upper curve) and ^ (lower 
curve) along with a numerical calculation of the 
quantity at two values of k (16 3 x32 lattice, (3 — 
5.9). It can be seen that for this quantity, not 
only is the unknown function of m not singular, 
it is approximately equal to 1. 

Putting the new action on a secure footing will 
ultimately require: 1) determination of the bare 
parameters of the action with mean field the- 
ory and full perturbation theory, 2) nonpertur- 
bative tests of the perturbative results, and 3) 
phenomenological tests of the resulting action in 
calculations of well understood physical quanti- 
ties. Not much of this program has yet been ac- 
complished. However, as argued above, at large 
values of ma, the new action (and even the Wil- 
son action suitably reinterpreted) can be viewed 



Figure 1. Normalization of the local vector cur- 
rent as a function of n. 



simply as unusual members of the general class of 
lattice actions proposed by Lepage and collabora- 
tors for NRQCD. Quite a bit is now known about 
the action for NRQCD. The discussion in points 
1) and 2) of Sec. 2.1 on operator classification is 
valid for any method for treating heavy quarks, 
including this one. The fact that the mean field 
corrections discussed in 4) reproduce the mass- 
dependent one loop corrections very well is en- 
couraging. 

Care will clearly be required in formulating nor- 
malization conditions which capture the most im- 
portant physics in both the relativistic and non- 
relativistic regions. (Identifying d£/dp 2 rather 
than M as the fundamental mass condition is ex- 
ample number one of these.) 



PHENOMENOLOGY 
AND T SYSTEMS 



OF THE J/ip 



Like all phenomenological lattice calculations, 
calculations of the properties of heavy quark sys- 
tems serve a variety of purposes. Quantities 
which are well understood experimentally, but 
which are very sensitive to lattice approxima- 
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tions are good tests of lattice methods (Sec. 3.3). 
Quantities for which the lattice approximations 
are well understood may be used to extract in- 
formation about the standard model (Sec. 3.1). 
A further purpose for lattice calculations is the 
delineation of the limits and the reasons for the 



successes of earlier models of hadrons (Sec. |4.3| ) . 

I will discuss calculations in the tfj and T sys- 
tems by Davies, Lepage, and Thackcr JL7| using 
NRQCD, and calculations in the ip system by 
the Fermilab group using Wilson fermions rein- 
terpreted as described in Sec. |2.2| p"8| - p2| ], and 
by UKQCD using Wilson fermions J23[| . (See also 
pi[.) Both groups studying the J/ip system used 
the 0(a) correction term 

Q 

5C = -ig-ipY.^F^ip (8) 

of Sheikholeslami and Wohlert ||. UKQCD 
used the tree level coefficient c = 1. The Fermi- 
lab group used a mean field improved coefficient 
c = 1.4 (see Appendix 

3.1. IS IP Splitting 

An excellent determination of the lattice spac- 
ing in physical units is provided by the spin av- 
eraged splitting between the lowest angular mo- 
mentum (/ = and 1 = 1) levels of the tp and 
T systems. (In the charm system, for example, 
M hc - (3M^ + M Vc )/4 = 458.6 ± 0.4 MeV.) The 
values of the lattice spacing obtained from this 
splitting do not differ dramatically from those ob- 
tained from other quantities, such as the p mass 
p(i[ or the string tension ^7|. It is the possi- 
bility of making improved uncertainty estimates 
that makes this an important way of determin- 
ing the lattice spacing. In quarkonia, error es- 
timates may often be made in several ways: by 
brute force (e. g., by repeating the calculation 
several lattice spacings), by phenomenological ar- 
guments, and by direct calculation of correction 
terms. Since determination of the lattice spacing 
is one of the key components of the determination 
of the strong coupling constant from low energy 
physics, it is important that these uncertainty es- 
timates be made rock solid. 

Preliminary results for this mass splitting were 
reported last year by the Fermilab group Jl8| , pj| 
and by Davies, Lepage, and Thacker IT?]]. This 



year, El-Khadra |2lj] reported further work done 
to check the corrections and error estimates given 
in Ref . jl8| . In jl8| , uncertainties due to an incor- 
rectly known quark mass and to 0(a) errors aris- 
ing from an imperfectly determined coefficient c 
in the 0(a) correction to the Lagrangian (Eq. ||) 
were taken to be less than 1% and omitted from 
the table of errors on the basis of the phenomeno- 
logical arguments. (The splitting is expected to 
be insensitive to small errors in the definition of 
the quark mass since it is almost identical in the 
ip and T systems. Likewise, in quark models the 
contribution of the a ■ B interaction, which domi- 
nates 8C nonrelativistically, to the spin averaged 
splitting is zero.) These arguments were checked 
this year by repeat calculation at several values 
of the parameters and found to be correct within 
statistical errors. 

The 0(a 2 ) errors were argued to be small be- 
cause repeat calculations at (3 = 5.7, 5.7, and 
6.1 yielded almost the same result for a(5 GeV) 
(see next section). An attempt was made to cor- 
rect for the small variation observed by extrap- 
olating to zero lattice spacing in a 2 . This ex- 
trapolation is not completely satisfactory, since 
the small a functional form is a messy combina- 
tion of 0(a 2 ) errors and perturbative logarithms. 
The way to improve this result, which has not 
yet be done, is to follow the example the NRQCD 
group |l7j and evaluate directly the contributions 
of the known correction operators to the splitting, 
thereby eventually obtaining zero measurable de- 
pendence on the lattice spacing. This group eval- 
uated the correction of the operators perturba- 
tivcly using the wave functions of the Richardson 
pj| ] potential model. For the T at (3 = 6.0, for ex- 
ample, they obtained the rather small corrections 
shown in Table [l| Wave functions directly calcu- 
lated by lattice gauge theory could also be used, 
eliminating the need for potential models. They 
are easy to calculate to high statistical accuracy. 
Fig. |^ shows the Coulomb gauge wave function 
of the J/ip meson calculated on a 24 4 lattice at 
j3 = 6.1. Statisitical errors are negligible at small 
separations. 

Similarly, the estimate of the finite volume cor- 
rection needs to be bolstered by calculating the 
meson Coulomb gauge wave functions on the lat- 
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the review of QCD in the 1992 particle data 
book. They are inconsistent with the most recent 
LEP determinations, which are around 0.120 and 
above. |§ 

The determination of a s from the 1S-1P split- 
ting currently consists of three separate elements: 
the determination of the lattice spacing, the de- 
termination of a physical coupling constant at a 
scale measured in lattice spacing units, and, for 
the time being, a correction for the absence of 
light quarks. As discussed in the previous sec- 
tion, the uncertainties in the determination of a s 
arising from the determination of the lattice spac- 
ing seem to be in good shape right now, and the 
path is clear to making them very solid. 



Figure 2. The wave function of the J/ijj meson. 



Term 


AM (IP- 15) 


% 


0(v 2 ) 


-11 MeV 


-2% 


0(a t ) 


13 MeV 


3% 


0(a%) 


-12 MeV 


-3% 




-24 MeV 


-5% 


Total 


-33 MeV 


-7% 



Table 1 

NRQCD corrections to the IP-IS splitting in the 
T system. Finite lattice spacing corrections are 
for P = 6.0. 



tice, and then calculating the overlap integral of 
the wave function with its periodically reflected 
image. 

3.2. Determination of a s from the IS— IP 
Splitting 

The most recent determinations of a s from 
the charmonium and bottomonium spectra us- 
ing NRQCD Jl7| and modified Wilson fermions 
p2| have error bars bracketing the region a s — 
0.103 - 0.114 GeV. They are somewhat below, 
but consistent with the world average given in 



3.2.1. Determination of the coupling con- 
stant. 

To determine the running coupling constant, 
one would like to combine the determination of 
the lattice spacing discussed above with a non- 
perturbative calculation of a physically defined 
coupling constant, for example defined from the 
static quark potential at a given, fixed momen- 
tum transfer like 5 or 10 GeV. 

Since the largest cutoff momenta for the exist- 
ing IP-IS splitting calculations was ir/a » 7.5 
GeV, it was not possible in the existing calcula- 
tions to determine the continuum limit of a phys- 
ical coupling defined at short distances. 

In Ref. |l8| a mean field improved perturbative 
relation, (Eq. [2?] in Appendix |b|) was used to ob- 
tain a renormalized coupling from the bare lattice 
coupling. This relation was tested over the past 
year on short distance Wilson loops. |3(J It did 
well, but not perfectly: the loops calculated by 
Monte Carlo were systematically a few per cent 
high. This suggests that Eq. 27 fixes most but 
not all of the pathological relation between the 
bare lattice coupling constant and physical cou- 
plings, and that it is better to obtain physical cou- 
plings from short distance quantities calculated 
nonperturbatively. Using short distance Wilson 
loops for this purpose (for example, via Eq. |28|) 
raises the values of the renormalized couplings by 
a few per cent over those reported in Ref. Jig]. 

It remains to be determined how much cou- 
plings defined from continuum quantities differ 
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from those defined from short distance quantities 
like the log of the plaquette. There is some rea- 
son to expect that this difference is small. The 
second order corrections to the short distance lat- 
tice static potential are within a few per cent of 
the continuum corrections. Likewise, Creutz 
ratios of Wilson loops up to six by six are quite 
well behaved when expanded with a coupling de- 
fined from Eq. |f.@ 

The Monte Carlo calculation of the static quark 
potential at a separation of one lattice spacing 
agrees to very high accuracy with perturbation 
using the coupling of Eq. ^8|. (See Sec. 4.1.) 
Therefore, a coupling constant obtained from the 
very short distance static potential will give re- 

A 
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suits almost identical to those using Eq. 
phenomenological method for estimating the con- 
tinuum coupling constant defined by the potential 
using short distance data has been proposed by 
Michael. H It has been used by UKQCD Q and 
by Bali and Schilling. It yields results which 
are quite close to those obtained with Eqs. |27] 
and |2^, and therefore to those which would be 
obtained directly from short distance lattice po- 
tential itself. 

A program to calculate explicitly a continuum 
coupling constant using finite size scaling has 
been proposed by Liischer et al.|Q To select the 
particular coupling to focus on, they propose the 
criteria that it: 1) be defined nonperturbatively, 
2) be calculable in perturbation theory, 3) be cal- 
culable in Monte Carlo simulations, and 4) have 
small, controllable lattice artifacts. These lead 
them to propose the response of the QCD vacuum 
to a constant background color-electric field to de- 
fine the coupling constant. The more phenomeno- 
logical choice of the static potential has poorer 
signal-to-noise properties in Monte Carlo calcula- 
tions approaching the continuum limit, and (they 
tell us) more difficult higher order perturbation 
theory. An SU(2) calculation has been completed, 
which yields results similar to those obtained with 
Eq. gjj. An SU(3) calculation is m progress. 

3.2.2. Correction for the effects of sea 
quarks. 

This is the greatest source of uncertainty in the 
results quoted above. This correction is the most 



phenomenological, and has the greatest likelihood 
of having a problem. Over the next few years, it 
will be removed by direct inclusion of the effects 
of sea quarks. 

The attempt to estimate the effect that the ab- 
sence of sea quarks has on this result is based on 
three assumptions. They are, in order of decreas- 
ing rigor: 

• When certain physics quantities are used to 
tune bare parameters in the quenched ap- 
proximation, the most important terms in 
the effective Lagrangian at the dominant 
energy scale for those quantities are given 
correctly. The effective action at other en- 
ergy scales including the scale of the lat- 
tice cutoff will be somewhat incorrect. In 
particular, if the effective coupling constant 
at the physics scale approximates that of 
the real world, the effective coupling at the 
short distance cutoff will be a bit small. 

• The most important term in the effective 
action for quarkonia is the static potential. 
The phenomenological success of potential 
models indicates that this assumption may 
be valid to around 25%. It is this assump- 
tion, which is certainly not valid for the 
light hadrons, that leads us to dare to try 
to make this correction for the charmonium 
system when we would not try it for the 
light hadrons. 

• The effects of light quarks on the static po- 
tential may be estimated by fitting char- 
monium data with a QCD based poten- 
tial model such as the Richardson potential 
once with the correct, rij = 3, [3 function in 
the potential and again with the quenched, 
rif — 0, (3 function. 

The final assumption is certainly a good one 
at short distances, which are responsible for most 
of the difference in the evolution of the coupling 
from the middle distance charmonium physics 
scale down to the lattice cutoff scale. It is 
also reasonable at the less relevant large distance 
scale, since the lattice quenched string tension 
and the string tension of Reggc phenomenology 
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are comparable. If, however, light quarks have 
a much greater effect on the potential in middle 
distances than they seem to at large and at small 
distpp ances, the assumption would fail. 

The naive expectation for the size of the cor- 
rection is 



Po ~ Po 

p7 =3 



20%. 



In Ref. 1 18 1, a perturbative calculation was used 
to bound the plausible size of the correction. This 
year the estimated correction was checked [^2| by 
fitting the charmonium spectrum with a poten- 
tial twice: once using a potential with the correct 
(3 function and once using a potential with the 
quenched (3 function. (See Sec. 1.2) The result 
was compatible with the one in Ref. This, 
however, is not so much an independent check of 
the previous estimate as another quantification of 
the assumption that sea quarks have no more dra- 
matic effects on the potential at middle distances 
than they seem to at large and small distances. 

3.2.3. Future prospects for determining a s 

Errors in the determination of the lattice spac- 
ing are already in good shape. 

The accuracy in the determination of the cou- 
pling constant needs further examination, but the 
calculations of Ref. suggests that the accu- 
racy to be expected of lattice perturbation theory 
is greater than that expected of QCD perturba- 
tion theory in hadronic phenomenology. In Ref. 
]30[ , discrepancies of about a 2 were typically ob- 
served in comparisons of first order perturbation 
theory with Monte Carlo calculations, and dis- 
crepancies of about a 3 in comparisons of second 
order perturbation theory. This amounts to only 
3-4% for calculations at the lattice cutoff at mod- 
erate /3's. In contrast, in QCD phenomenology, 
an accuracy of 10% is often taken to be opti- 
mistic. One difference may be that the lattice 
calculations of most interest are often quadrati- 
cally divergent integrals dominated by momenta 
of the order of the relatively well-defined lattice 
cutoff. They thus differ from calculations of un- 
ruly hadrons in collision, which insist on interact- 
ing on a wide range of momentum scales, piling up 
large logarithms from a nasty variety of sources. 



The aspect of the current determination of a s 
which makes it no better than any other existing 
determination is the use of potential model argu- 
ments to estimate the effects of the absence of sea 
quarks. This is quite analogous to, for example, 
the phenomenological treatment of higher twist 
and fragmentation effects in determinations of a s 
in deep inelastic scattering. All of the existing 
determinations have some phenomenological as- 
sumptions built into them. The difference is that 
the potential model estimate of quenched correc- 
tions will certainly be eliminated by brute com- 
puter force (if not by the use of more intelligent 
methods) over the next few years, resulting in a 
determination far more accurate than any of the 
existing ones. 

3.3. Hyperfine Splitting and Leptonic 
Width 

These two quantities are very straightforward 
to calculate on the lattice, and are good phe- 
nomenological tests of how well we understand 
the parameters of the quark action. The hyper- 
fine splitting Am( J/ip — rj c ) and leptonic decay 



amplitude Vw, 



)J have been calculated in 



the J/V> system by the Fermilab group and by 
UKQCD. They have been calculated in the T sys- 
tem by Davies, Lepage, and Thacker.[jl7j 

The potential model formula for the hyperfine 
splitting is 



Am(J/^ - Vc ) = 9 ^ l 2 > (0)| 2 



(9) 



It arises from a coupling of the spins of the quarks 
to transverse gluons. It therefore should be ex- 
tremely sensitive to the value of the correction 
coefficient c. (It is not clear a priori whether to 
expect strong sensitivity to the quark mass, since 
^(O)! 2 should rise with the quark mass.) 
The leptonic width to leading order is 

mi 

Nonrclativistically, the leptonic decay amplitude 
is therefore simply the wave function at the origin, 
^(0), properly normalized. This quantity should 
be quite sensitive to the mass of the quark. 

Before comparing existing lattice results with 
experiment, we need to estimate the accuracy 
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to expect in the quenched approximation. Both 
quantities are proportional to |^(0)| 2 , the prob- 
ability for the quarks to be at the same point 
(within one Compton wave length, say) and so 
are obviously short distance quantities. With lat- 
tice parameters tuned to obtain the correct IP-IS 
splitting, the coupling constant and potential at 
short distances will be too weak. An analysis like 
the one referred to in Sec. 3.2.2 yields 



(0)/ 1 

as '{m c ) 

(3)/ \ 

Q-s {m c ) 



0.81 ±0.06 



(11) 



The incorrect weakness of the quenched po- 
tential at short distance may also yield a weak- 
ened wave function at the origin. El-Khadra has 
checked for this effect 21 2^] using the Richard- 
son potential 



V(q 2 ) = C f 



where /?q ^ = 11 — 2n//3 . Fitting the charmo- 
nium spectrum with once with rif = 3 in the (3 
function parameter, and again with rif = 0, she 
found that, indeed, The ratio of the wave func- 
tions at the origin was 



4-7T 



1 



q 2 ln(l + 9 2 /A 2 ) 



(12) 



fr(°)(0) 

*(3)(0) 



0.86 



(13) 



This reduces our expectation for the hyperfmc 
splitting from the experimental result Am(J/fi — 
,y c ) cx P = H7.3 MeV to around 



Am( J/tp - T]c ) ( i ucnchcd ~ 70 MeV 



(14) 



These considerations reduce our expectations for 
leptonic matrix element by a smaller amount, 
from the experimental result V^ xp = 0.509 
GeV 3 / 2 to 



•4> 



inched 



0.438 GeV 3/2 



(15) 



The results of UKQCD and Fermilab for the 
hyperfmc splitting are shown in Fig. § UKQCD 
set the value of the quark mass to obtain the ex- 
pected energy eigenvalue in the transfer matrix. 
In light of the arguments on the interpretation of 
Wilson fermions at large quark masses in Sec. 2.2 



Figure 3. Ara(J/i/> - rj c ) calculated by UKQCD 
(diamond) and Fermilab (squares) compared with 
the physical result (upper star) and an estimate 
of the quenched corrected result (lower star). 



the Fermilab group took this as unreliable and at- 
tempted to fix the quark mass by demanding that 
the leptonic width be correct. The UKQCD result 
is slightly below the result expected on the ba- 
sis of the quenched correction. This is consistent 
with the mean field expectation that quantum 
corrections boost the required value for the co- 
efficient of the correction term. Their results are, 
however, much closer to the physical answer than 
earlier Wilson fermion calculations with no 0(a) 
correction. [^5| The Fermilab results are slightly 
above the quenched expectation, but perhaps not 
very significantly in light of the uncertainties in 
the quenched correction and the statistical errors. 

4. THE STATIC QUARK POTENTIAL 

High accuracy results for the static quark 
potential were reported this year by Bali and 
Schilling @|3 and by UKQCD §§. Fi S- 
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Figure 4. The heavy quark potential, calculated 
on the lattice in the quenched approximation. 



|4] shows the potential calculated by Bali and 
Schilling in the quenched approximation on a 32 4 
lattice at (3 = 6.4. The solid line is the fit to a 
Coulomb plus linear potential 

V(R) = V -0.277(28)/# + 0.0151 (5)R, (16) 

which fits quite well for R > 2\/2. Comparison 
of results on 16 4 with results on 32 4 lattice indi- 
cated that finite volume results are small. A plot 
of results from f3 = 6.0, 6.2, and 6.4 with physi- 
cal units set by the string tension indicates good 
scaling behavior. 

4.1. Short distance behavior. 

At such a large (3 we should expect the short 
distance part of the potential to agree very well 
with perturbation theory, and this is the case. 
I checked the value of V(l) given in Ref. |2^] 
against perturbative results for 32 4 lattices sup- 
plied by Urs Heller 37 1. Using the "measured" 
coupling constant defined by Eq. |H, perturba- 
tion theory agreed with the Monte Carlo data to 
within about 1%, perhaps fortuitously accurate, 
but still impressive. (This incidentally illustrates 
that the potential is a natural candidate on the 
lattice as well as in the continuum to define im- 
proved coupling constants. The coupling of Eq. 
|2§| was suggested mostly because the plaquette is 
easy to measure and universally available.) 



Since perturbation theory agrees so well with 
the Monte Carlo calculation of the potential, and 
since perturbation theory implies a coupling con- 
stant rising with increasing R, it would be inter- 
esting to attempt to fit the data with an asymp- 
totically free Coulomb plus linear potential. The 
size of the fit Coulomb term (0.277(28) / R) is 
quite close to the subleading long distance behav- 
ior of the potential (tt/(12R) = 0.262/i?). How- 
ever these two similarly-sized effects have nothing 
to do with each other, and we are not guaran- 
teed that, for example, the perturbative Coulomb 
term does not rise above 0.28 before the poten- 
tial settles back down to its asymptotic form. A 
fit with an asymptotically free Coulomb plus lin- 
ear potential, for example a modified Richardson 
potential |58|, might help to start exploring the 
extent to which the data support or rule out such 
speculation. 

It is easy to convince yourself with a ruler that 
values of the string tension obtained by the fit are 
completely plausible. However, if the changeover 
from the perturbative Coulomb potential to the 
nonperturbative long distance 1/R term is more 
complicated than we hope, a larger than expected 
middle distance Coulomb term could be contam- 
inating the obtained string tensions more than is 
obvious from the current analysis. 

4.2. Long distance behavior. 

The good scaling of the potential when the 
physical scale is set by the string tension has al- 
ready been mentioned. Good asymptotic scaling 
of the string tension in terms of a physical cou- 
pling such as Aj^g is also observed in the new 
data. (Good means to perhaps 20%.) Folklore 
to the contrary was based on the the search for 
scaling in terms of the bare lattice coupling con- 
stant. The bare coupling has a highly patholog- 
ical, but reasonably well-understood relationship 
to well-behaved physical coupling constants. It 
was pointed out long ago |3j],|t| that decent scal- 
ing is observed in terms of an effective coupling 
constant defined from the plaquette. It was em- 
phasized in Ref. |?(| that such coupling con- 
stants are simply very close relations of the fa- 
miliar physical coupling constants such as ay and 
a Ws °f perturbative QCD. 
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Figure 5. V^/^Ws as a function of the lattice 
spacing. The upper curve was obtained from the 
bare coupling constant. The lower curve was ob- 
tained from an effective coupling constant. 



Fig. ^ (from Ref. j4(]]) shows new and old data 
for yfa/ k-jj-g plotted as a function of aAj^g. The 
upper curve was obtained via the bare coupling 
constant. The much better behaved lower curve 
was obtained via the effective coupling 

3(l-|rr(LQ) 



a eff 



(17) 



47T 

(Ref. [[30| advocates Eq. ^8|, the logarithm of 
Tr(U), for this purpose on the grounds that log- 
arithms of Wilson loops have better perturba- 
tive behavior than the loops themselves.) Only 
about 20% deviation from asymptotic scaling is 
observed over the range of the data. Part of that 
deviation is certainly perturbative, since the use 
of another reasonable perturbative scheme, Eq. 
g|[ changes the amount of deviation to 10%.[p0| 
The fact that the ratio of the deconfinement tem- 
perature to the square root of the string ten- 
sion scales better than ^fa / 'A^g-^ is another in- 
dication that the deviation is more likely to be 
connected to the determination of the coupling 



Figure 6. The lattice quenched heavy quark po- 
tential (top curve) and the potentials of the Cor- 
nell model and the Richardson model. 



constant than to the string tension. Because the 
short distance behavior is a mixture of pertur- 
bative logarithms and 0(a 2 ) errors, the extrapo- 
lation in a is not completely satisfactory and it 
is important to sort the origin(s) of the discrep- 
ancy: perturbation theory, 0(a 2 ) errors, or mea- 
surement errors. However, the downward trend 
seems clear and the estimate 



V / °/ A TTs = 1-75 ±15% 
seems reasonable. 



(18) 



4.3. Comparison with potential models. 

One useful task of first-principles calculations 
is to support or destroy earlier phcnomcnologies. 
For example, it would be nice to be able to un- 
derstand if there is a reason that nonrelativistic 
quark models for the light hadrons work unrea- 
sonably well. It is more straightforward to put the 
success of potential models of heavy quark sys- 
tems on a rigorous footing using lattice methods. 
These systems are nonrelativistic and it is not 
surprising that a nonrelativistic treatment yields 
rather accurate results. 

In Fig. H the potential obtained by Bali 
and Schilling is compared with the potentials of 
Eichten et al. M] and Richarson in the region 
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0.1 fm < R < 1.0 fm. The string tensions of the 
lattice and the phenomenological potentials are 
similar, but the Coulomb term required by phe- 
nomenology is about 1.8 times as large as that 
yielded by the quenched lattice, seemingly a large 
discrepancy. The phenomenological potential is 
very well known in this region between 0.1 and 1.0 
fm. Fits to the spectra of the charmonium and 
bottomonium systems with a wide variety of plau- 
sible and implausible functional forms yield po- 
tentials which differ by only a few per cent in this 
region. On the other hand, the quenched lattice 
potentials are also rather convincing, especially 
at short distance, so what accounts for the dif- 
ference? First, we expect the quenched Coulomb 
coupling to be a bit smaller than the true QCD 
coupling constant at short distances because of 
the incorrect (3 function of the quenched approx- 
imation. (See Sec. 3.2.2.) This effect is in the 



right direction and is expected to be of order 



n + — 3 



20%. 



Second, the phenomenological potentials clearly 
parameterize some of the effects of higher order 
relativistic corrections. These are roughly ex- 
pected to be of order v 2 /c 2 ~ 25% for charmo- 
nium. Some of these clearly have the effect of 
strengthening the attraction of the quarks, but 
a complete analysis of the spin-independent rel- 
ativistic corrections in potential models does not 
exist. 0] A combination of these two effects could 
thus easily explain as much as 1.5 out of the dis- 
crepancy of 1.8. A preliminary conclusion: there 
is an interesting puzzle in this discrepancy, but 
no cause for alarm. 

5. SUMMARY 
Static Potential. 

• At short distances, the potential agrees with 
perturbation theory to a few per cent. 

• The string tension exhibits two loop asymp- 
totic scaling to an accuracy of 20%. 
^/a/Aj^g is in the range 1.55-1.95. 

ip and T systems. 



• The hyperfme splitting and leptonic widths 
provide good phenomenological tests of lat- 
tice methods. 

• The spin averaged IP- IS splitting provides 
a very good determination of the lattice 
spacing in physical units. Combined with 
a lattice determination of the renormalizcd 
coupling, it gives a determination of the 
strong coupling constant which at present 
is of comparable accuracy to that of con- 
ventional determinations. When the effects 
of sea quarks are properly included, its ac- 
curacy will be much better than any current 
determination. 

Technical developments. 

• Lattice perturbation theory works very well 
when renormalized coupling constants are 
used. 

• Minor changes to the actions of Wilson and 
of NRQCD may make it possible do cal- 
culations with a unified formalism at any 
value of the quark mass, as long as the 
three momentum is small. This will imply a 
reinterpretation of calculations with Wilson 
fcrmions at large quark mass. 
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A. NOTATION 

We use the forward, backward, and symmetric 
finite difference operators 



A-lP n = tp n - C^-^Vn-Ai 
A^Vn = ^ ~t -V'n- 



(19) 
(20) 

(21) 
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The analogous continuum covariant derivative is 
denoted 



(22) 



The standard relation between the bare mass 
mo and the hopping parameter k is 



to 



-i-4. 
2k 



(23) 



When considering nonrelativistic fermions, we 
decompose the four-component Dirac field as two 
two-component fields 



x 



(24) 



and take as our representation of the Euclidean 
gamma matrices 



To 



1 

-1 



7i 



<J t 
<Ti 



(25) 



B. RESULTS FROM LATTICE PER- 
TURBATION THEORY 

This section summarizes results from Ref. [ fjo) 
which have been used in the text. 

B.l. A sequence of improved coupling con- 
stants 

In Ref. |H (1990) it was argued that lattice 
perturbation series are much more convergent and 
agree better with Monte Carlo data if they are 
expressed in terms of a physical running coupling 
evaluated at a carefully chosen scale. A good one 
is ay, the one defined by the static quark poten- 
tial: 



= — +f3 ln(— ) - 4.702. 

OtVW Otlat a <l 



(26) 



The arguments were analogous to those leading 
from the M S to the MS scheme in dimensionally 
regularized QCD. 

In Ref. [0] (1991) it was noted that the bulk of 
the correction coefficient in the previous equation 
is accounted for by a simple mean field argument. 
The coupling constant is enhanced at one loop 
by a coupling to the expectation value of the pla- 
quette induced by the higher order terms in the 
Wilson action. Higher order analogues of this one 



loop effect certainly exist. This suggests that an 
effective coupling constant which incorporates a 
Monte Carlo calculation of the plaquette expec- 
tation value, such as 



1 



a v (q) 



(TrU ) 

\ SL + l n (— )- 0.513, (27) 



OLlat 



7T 

aq' 



may yield improved accuracy. 

Over the past year (1992) we have tested this 
assumption by calculating a variety of short dis- 
tance quantities using the mean field improved 
coupling constant and by Monte Carlo. ^| We 
found that, while using Eq. ^ significantly im- 
proved agreement between perturbation theory 
and Monte Carlo, the Monte Carlo results tended 
to be systematically slightly higher then the per- 
turbative results. This suggests that a coupling 
defined directly from any of the Monte Carlo cal- 
culated quantities would yield improved predic- 
tions for the others. A particularly simple one is 
the coupling defined from the log of the plaquette: 



1 



a v (3Al/a) 



1 

a eff 



1.19 

-in 



31n(±Tr([/)) 



1.19 



(28) 



(The scale of the running coupling arises from an 
estimate of the typical momenta of gluons in the 
calculation of the logarithm of the plaquette. [po[) 

B.2. Mean field improvement of operators 

The mean field argument leading to Eq. ^ may 
be summarized as follows. The naive classical 
relation between the lattice and continuum gauge 
fields 



I7 M (ar) = e iagA ^ -> 1 + iagA^x). 



(29) 



is spoiled by tadpoles arising from the exponential 
form of the lattice representation of the gauge 
fields. Quantum fluctuations do not lead to an 
average link field close to 1.00000 as implied by 
Eq. E9, but to something more like 



U^(x) — > u Q (1 + iagA^x)), 



(30) 



where uq, a number less than one, represents the 
mean value of the link. In a smooth gauge, the 
Monte Carlo link expectation value can be used 
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as an estimate of uq. A simple, gauge-invariant 
definition is 



w 



plaq 



.1/4 



(31) 



Other definitions based on re c or the static quark 
self-energy may be used to fine tune mean field 
predictions in particular situations. 

If naive definitions such as Eq. ^9] are used 
to relate lattice and continuum operators, large 
corrections will appear in quantum corrections. 
Much better behavior of loop corrections is ob- 
tained by taking U^(x)/uq as the lattice approxi- 
mation to the continuum field. This implies that 
the lattice action 



=duon 



„2„4 



Tr([/ plaq + h.c). 



(32) 



will approximate closely the desired continuum 
behavior. This is the usual lattice action if we 
identify 



~2 



gfat/uo = gfat/(kTr(U p i aq )). 



(33) 



The perturbative result, Eq. |27], explicitly verifies 
that g 2 is a closer approximation to a standard 
continuum expansion parameter than gf at is. Gil 
The same considerations lead to the result that 



re = KUq 



(34) 



produces a more continuum-like bare mass [fh — 
— 4) and smaller quantum corrections in oper- 
ator renormalizations than does re. 

Mean field arguments may be used both to es- 
timate perturbative predictions when the pertur- 
bative predictions are unknown, and also to im- 
prove known predictions. Just as we did with 
the coupling constant in Eq. |27], we can improve 
perturbative predictions for operators involving 
quark fields by substituting the Monte Carlo cal- 
culation of u in Eq. |34| and including in the 
perturbative prediction only that part remaining 
after the absorption of uq into re. A possible fine 
tuning in this case is to obtain uo from re c rather 
than from the plaquette. This was the procedure 
used in obtaining Eq. [?]. 

The cloverleaf approximations to F^ v used in 
the Wohlert-Shcikholcslami 0(a) correction to 
the Wilson action |25| and in the magnetic spin 



coupling of NRQCD |Ll| contain four links each. 
The quark wave function normalization contains 
one link. We therefore expect the naive coeffi- 
cient of this operator to undergo quantum cor- 
rections of roughly a factor of Uq 3 , or about 
1 + 0.25g 2 if we use the plaquette to estimate 
Uq. This is in agreement with an unpublished 
thesis calculation of Wohlert, 1 + 0.27g 2 .J|| Us- 
ing the plaquette calculated by Monte Carlo to 
estimate the correction term yields a factor of 
(±Tr([/ plaq ))- 3 / 4 - 1.4 - 1.5 for (3 around 6.0. 
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